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eleaent  of  surface  area 
drag  coefficient,  Eq.  3.4 
lift  coefficient,  Sq.  3.5 

partial  drag  ana  lift  coefficients  proportioned  to  term 
indicated  in  subscript.,  Sqe.  3.4,  3.5 

energy  fluxes  incident,  reflected  and  diffusely  reflected 
from  dA  , Eq.  2.7 

drag  force  per  unit  length  of  cylinder 
lift  force  per  unit  length  of  cylinder 
force  on  cA 

component  of  di?  due  to  incident  molecules 
component  of  dh  due  to  reflected  molecules 
molecular  velocity  distribution  function,  Eq.  2.5 
Maxwellian  distribution,  Eq.  2,4 

modified  Sessel  functions  of  sero  and  first  order 
thermal  conductivity 
Mach  number 
molecular  mass 

number  of  molecules  incident  on  <dA  per  unit  time 
pressure 

pressure  on  dA  due  to  incident,  reflected,  and  diffusely 
reflected  molecules,  respectively,  Eq.  2. 1C 

normal  stress  deviation  terms,  Eqs.  2.5»  2.6 

heat  flux  terns,  EqB.  2,5,  2.6 

cylinder  radius 

gas  constant 

molecular  speed  ratio,  Eq.  3.3 

molecular  speed  ratio  refered  to  cylinder  temperatures,  Eq,  3.3 
( Sy^  assumed  equal  to  S in  cor,  puts  ties  a) 

gas  tenperature 

temperature  of  incident  molecules,  Eq.  2,13 
cylinder  tenperature 

molecular  velocity  component  in  X , y ,Z  direction,  Fig.  1 
gas  velocity  components 

molecular  velocity  components  in  X' , y 1 , 2.1  directions,  Eig.  1 
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coordinates  rtvfered  to  flow  direction.  Fig.  1 
coordinates  refered  to  surface  alsoent,  Fig.  1 
thersal  aceconodation  coefficient,  5q«  2.7 
ratio  of  specific  heats 
viscosity  coefficient 

angle  between  ourfacs  noraal  and  flow  direction,  Fig.  1 
gas  density 

reflection  coefficients  for  tangential  and  noraal  scwOatus 
tmnsfer,  Xqa.  2.9,  2*10 

shear  stress  due  to  incident,  reflected  and  diffusely  re- 
flected aolecules,  respectively,  Fq.  2.9 

shear  stress  terns,  Eqa.  2.5,  2,6 
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1.0  ireROPtCIIQN 

The  forces  exerted  open  a body  in  fr^e  molecule  floe,  i.e„,  on  a 
body  whose  characteristic  dimension  is  small  compared  to  th>3  molecular  mean 
free  path,  has  been  determined  theoretically  by  Teles  (Ref,  1) , Stalaer  and 
XuricJr  (Ref.  2)  and  others  (Refs.  10  to  12).  Experiments  by  Stalder,  Good- 
win and  Creager  (Raf.  3),  and  by  Sans  and  Eeterwann  (Ref.  4)  iwu'o  provided 
general  confirmation  of  the  theoretical  developments.. 

These  invest! gat Iona  have  been  confined  to  the  case  is  which  the 
gas  flow  past  the  body  ia  in  laxwellian  equilibria*.  This  restriction  la 
satisfactory  for  free  flight  considerations  but  not  for  what  way,  perhaps, 
prove  to  be  oca  of  the  most  important  applications  of  free  molecule  flow 
analyses.  This  ia  wind  tunnel  weesurement  of  the  velocity,  and  perhaps 
ether  flow  variables,  of  a gas  stream  by  weans  of  small  "free  molecule  flew* 
probes.  The  possibilities  are  especially  promising  in  connection  with  hyper 
sonic  or  supersonic  wind  tunnels  operating  st  low  test  section  densities. 
These  probee  mi^it  be  of  either  the  beat  transfer  or  ssro^ynumic  force  type, 
or  perhaps  both.  It  is  clear,  however,  that  for  probing  boundary  layers  or 
shock  waves  or  other  characteristic  parts  of  a flow  field,  the  instrment 
will  often  be  exposed  to  a gas  stream  which  is  not  in  Maxwellian  equilibritmi 
The  occurence  of  a large  shear  stress  in  a boundary  layer,  for  example,  pro- 
duces a departure  of  the  molecular  velocity  distribution  from  the  Maxwellian 
As  will  be  seen  below,  this  non -uniformity  gives  rise  to  additional  force 
terms- 


Approximate  analyses  of  the  forces  on  bodies  in  such  non~“hx- 
welllan  free  molecule  flows  have  been  made  by  Einstein  (Ref.  0)  and  Sp- 
stein  (Ref.  6)  in  connection  with  the  radiometer  problem.  In  their  dis- 
cussions the  non-uniformity  was  confined  to  the  case  of  a local  heat  flux, 
and  in  any  case  the  results  were  only  of  order-of-magnitude  velidity. 

In  the  present  report  a wore  detailed  anelyais  of  the  free 
molecular  flow  of  a non-uniform  gas  is  made.  The  results  are  general 
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enough  to  include  arbitrary  (sjmII)  streos  and  fasot  flux.  The  body  geometiy 
has  been  confined  to  that  of  a right  circular  cylinder  oriented  perpendicu- 
lar to  the  gas  flow  velocity,  and  it  la  assumed!  the  beef's  thermal  conduct- 
ivity ie  so  large  that  it  ie  at  uniform  temperature.  Only  the  aerodynamic 
force  characteristics  are  determined.  The  eAtensionc  to  the  determination 
of  heat  transfer  characteristics,  to  other  body  geometries  and  to  the  ease 
of  very  low  internal  conductivity,  will  be  presented  in  subsequent  reports. 
The  aethod,  however,  will  be  apparent  frea  tbs  present  peper. 


2.0 


2.1  Xk?..Ia  si  6 vP  t.  reel  u 

The  general  Method  followed  is  a direct  extension  of  the 

Method  for  the  Maxwellian  case  used  by  Tsiec  (Ref.  1)  and  Stalder 

and  Zurich  (Hef.  2),  and  waa  suggested  by  Goodwin.  The  force  per 

y . 

unit  trs*  dy  dA  on  a surface  element  dA  ia  broken  up  into  a 
part  d ’t~ j /d A due  to  the  incident  molecules,  and  a part  d^^/d A due 
to  the  reflected  or  re-emitted  moleculee.  Total  force  characteristics 
aro  obtained  by  integration  over  tba  surface  of  the  body. 

The  foree  dFj  y&A  is  analysed  in  terms  of  the  component 
p;  normal  to  the  surface  (pressure)  and  the  component  ?,'■  tangential 
to  the  surface  (shear  stress).  See  Fig.  1.  These  quantities  are 
obtained  in  teres  of  the  solecular  velocity  distribution  function, 
f (q,v,  W)  , which  is  the  n’.mber  density  per  unit  volwe  of  phase 
space  of  molecules  with  velocity  components  U,V,W.  One  has, 

*°  Z 

u‘  f d u d v' d w' 


u‘  v‘  f dud  vldw‘ 


(2.1) 


(2.2) 


where  fn  is  the  molecular  mass.  The  velocity  components  U , V, W and 
u'iV'jW'aro  related  by  (see  Fig.  1), 

u = - u‘  CCS  0 + v‘  sin  6 

v «~u'  sin  9 - v'  cos  8 > 

w * w‘ 


(2.3) 
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where  6 ic  the  angls  betweer  the  (outward;  norcsl  ti  the  surface  and 
the  X (flow)  direction. 

The  function  f 1b  given  in  the  Maxwellian  <saee  by 
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f (u,v,w)  -f  “=  nn'(2TT  pt)  ^ e ^ 

whereto,  ',U?R  are  the  ga£  density,  temperature,  flow  velocity  and 
the  gas  constant,  respectively.  For  the  non -Moxwe Ilian  case,  one  usee 
the  distribution  function  considered  by  Maxwell  (Ref.  7)  and  Grad 
(Ref . 8)  b 

i+  rpfrr!  Py //  + Pzz^2 


f = T 


■ 1 7rxy  (y-U)  V -2  r z (u  -U)  \N  -Z~yz_  \jsN 
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where  px><  9 ><  * rtc’  ar®  7l5coue  stress  ai  d heat  flux  com- 

ponents. These  quantities,  together  with yO?  U and  "T , are  treated  as 
constants  throughout  the  region  occupied  by  the  cylinder.  For  the 
fcquilihrirm  case,  and  for  stall  gradients  of  velocity  and  temperature, 
xhase  tents  are  related  to  the  velocity  and  temperature  by  the  Haviar- 
Stokee  relations 
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ror  rap  id  processes  ve.g.,  shock  waves)  or  for  large  gradients  of 
tenperature  and  velocity  more  complicated  relations  exist;  see  Ref. 
S for  details.  In  the  present  report,  the  results  will  be  obtained 
in  terms  of  the  variables  listed  in  Eq,  2,5,  so  thnt  they  arc  indc - 
pendent  of  whether  or  not  Fqs,  2,6  are  valid.  It  should  also  be 


(2.4) 


(2.5) 


(2,6) 


nentioned  that  for  sufficiently  non -uniform  situations  the  approximation 
of  Sq.  2.5  itself  becomes  poor.  Additional  teres  in  the  expression  for 
the  molecular  velocity  distribution  would  ther.  be  required  and  would  con- 
tribute to  the  gross  force  exerted  upon  the  body.  Ths  relative  depen- 
dence upon  the  parameters  b^C*  * howr>Hi'>  "O'uid  bo  unchanged. 

.2  Iks  faglgatofljtelsgula 

~Tjt 

The  determination  of  di^  . the  force  on  d A due  to  the  re- 
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v.ire*  a soeoif ication  of  the  nature  of  the  intor- 
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action  of  the  impinging  aolecules  with  the  surface.  The  trwtsont  of 
this  question  will  be  somewhat  generalised  in  thia  roport  free  the  usual 
treatment  of  maxwell  (Kef.  7}  which  has  bean  followed  in  the  subsequent 
sorb  of  Tsien  (Raf.  1),  Stsldsr  and  2uridt  (Raf,  2) } and  others.  Tra- 
ditionally, the  interaction  has  bean  described  in  tares  of  two  para- 
meters, the  thermal  accomodation  coefficient  c*  and  the  specular  re- 
flection coefficient  (T  . The  quantity oL  ia  defined  by: 

oc  = d^'  ~d^r 

dt,,  -dE* 

where  dE  ( , dEj.arddE^are  the  energy  fluxes,  respectively,  incident 
on  and  remitted  from  the  surface,  and  the  flu-  which  would  be  re- 
mitted if  all  incident  molecules  were  remitted  with  a lares Ilian 
distribution  corresponding  to  the  surface  ves  pc  miners u quan- 

biby  oL  then  is  a measure  of  the  degree  to  which  the  incidert  mole- 
cules are  =acct*ouatwd*  to  the  surface.  Measured  values  of  oL  are 
usually  slightly  less  than  unity  ('Ref.  9) . It  will  be  observed  thsi 
all  quantities  cn  the  right  side  oi‘  7q.  2.7  are  well-defined  and  can 
bs  either  calculated  or  ■esBurc-d. 


The  definition  of  (T  , however,  is  somewhat  less  satisfactory. 
It  is  supposed  that  the  incident  wcleeules  are  either  reflected  “spec- 
ularly* (i.e. , with  simple  reversal  of  normal  velocity)  or  :<iiff■uEe:^r", 
(which  for  eo»6  applications  can  be  just  "randomly*,  but  eventually 
must  be  defined  more  precisely  as  Maxwellian  corresponding  to  a toa- 
peiatureT^.  which  ie  not  necessarily  the  Base  anT^j).  The  coeffi- 
cient 0 is  then  defined  as  the  fraction  of  diffusely  reflected  aole- 
cul.se,  I - (J  being  the  fraction  which  are  specularly  reflected.  The 
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definition  in  this  fora  of  the  quantity  0 thus  becoees  sear,  in  glee  a 
unless  the  actual  interaction  with  the  surface  is  of  the  type  speci- 
fied, Ab  measured  valwe  of  (T  have  actually  been  determined,  how- 
ever, this  difficulty  hgg  been  avoided  because  only  groee  tangential 
forces  are  dealt  with.  Thus,  according  to  the  definition  of  (T  , one 
has  g 


6~  ~ • ^ r 

^ Tj 

where  £•  and  are  the  incident  and  reflected  fluxes  of  tangential 
aooentua  (shear  stress).  For  ail  practice!  purposes,  Eq.  2.8,  rather 
than  the  specular-diffusa  aodel , has  been  the  defining  equation  for 
Q . The  two  definiiiouu  are  in  agrees ant  if  the  interaction  aetu- 
ally  la  of  the  specular-diffuse  type  (which  it,  of  course,  almost 
certainly  is  not),  otherwise  Fq,  2,8  defines  (T\  Eq,  2.8  nay  be 
written  slao,  for  ayraaatry  with  Eq.  2,7,  In  the  fors: 


(2.8) 


(J-  as  ^ t - C r (2,9) 

/7~'.  — '7^ 

U \ L vv 

where  'f ^ is  defined  as  the  tangential  moaentm  flux  which  would  be 
reflected  from  the  surface  if  all  incident  oolecvJ.es  were  re  emitted 
■ith  a Haxwellian  distribution  corresponding  to  the  surface  tempera- 
ture"^. Clearly  0 , 


Thus  the  difficulty  inherent  in  the  original  definition 
of  O’  is  avoided,  insofar  as  the  determination  of  tangential  forces. 
However,  it  occurs  again  in  the  determination  of  aurukil  forces,  and 
again  recourse  has  been  taken  in  Refs.  1 and  2,  for  example,  to  the 
specular-diffuse  aodel.  If  thy  actual  interaction  is  different, 
then  the  determination  of  heat  transfer,  tangential  force  and  normal 
force  in  t ena  of  .lust  the  two  par-voters  o*-  and  0~  would  yield  results 
which  one  would  expect  to  be  mutually  inconsistent.  This  difficulty  can 
be  avoided  if  an  additional  parameter  U ; say,  is  defined  or 

- Pr 


(J  = 


r\ . 

K I 


Pi  -pw 


(2.10) 


where  p.  , pr  andp^are  respectively  the  incident  and  reflected 
fluxes  of  normal  r.owentum  (pressure)  and  the  flux  which  would  bo 
reeeitted  if  all  the  incident  nolacules  sere  reeaitted  with  a Max- 
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weliian  distribution  corresponding  to  the  surface  temperature""^,., . 
The  quantity  p is  given  by  (Ref.  2) 


o — \N  . nn  \i 
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(2.11) 


where  is  the  masber  of  molecules  per  unit  tiae  incident  on  a unit 
surface  area.  It  is  givsn  to 

,oo  O«0 

/ / X 


Nw  _ / / f u'f  d U1  d v'  d w' 

i L L 
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(2=12) 


Fy.*  completely  specular  reflection  one  has  od.  — (j  — (5~  O, 
For  completely  diffuse  reflection,  w = 0”**  1 . For  the  hypo- 

thetical partly  diffuse,  pertly  specular  interaction,  and  f„?  a sur- 
face at  rest  with  respect  to  the  gas,  one  would  have 


(J 
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> oc 


(2.13) 


However,  in  general,  for  actual  physical  interaction  the  three  pero- 
netera  would  be  independent . They  arc  sufficient,  however,  to  deter- 
mine free  molecule  flow  heat  transfer  and  aerodynamic  force  character- 
istics . 


2 . 3 The  C cabined  Expressions 

The  total  force  on  a unit  of  surface  area  is  analysed  into 
a normal  component-  p and  a tangential  component  'f  given  hf 

P * P',  + Pr 

^ ~ _ *>> 

h - tj  lr  J 

With  the  help  of  Eqs.  2.9  and  2.10  these  reduce  to 
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P - (2-  cr1)  p.  + P>, 
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r = C*  r: 


(2.15) 
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Tn  of  \ . xjb ing  Bqa,  2.11,  2.17.,  2.1  and  2.2,  one  has 

OO  <H£~> 


^ rn  / / /^f  Hz-G'Ju  ^u]d,dv!dw' 

-Sd  <T'  > 


u‘  v1  *f  d u d v1  a1  \aj 


The  total  force  on  e ■nr.it  length  of  the  cylinder  is  obtained 
in  teroa  cf  the  drag  force  r„  and  the  lift  force  F,  , These  are  given 
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p ^ j'  (-psine-f  cos©)rde  j 

Combining  Bqs.  2.17,  2.16,  2.5  and  2.3,  oce  has  finally, 
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5 RT 
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QO  O 
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e Z RT 


rW(-u'co$©  + v'sin 0-U)  + R,y(-u's,in 0 -v'coe 9)  ^Pzz^  (2.19) 

-?tr  (-u'cosy  Wsi  n 9 -U)(-u'  sin  & W COS  6)  -Z  1“ , (-u'cose  +v!s'in  6-Ly  vj 
" xy v \ 
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,y~  ■ \ 5 RT  j 

Js'x 

(q  _(-u,cosetv'&ine-^+qy(-u,9ir.9-v,cose)+qti«')  li  rdu'dv'dw'de 

T Tj 


Performing  the  integration  indicated  in  Bqa.  2.13  and  2.19,  one 
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(3.2) 


/,  JST  fSS]  + Mi-lt  (Sl)  + T (W 

ebere  I and  X are  the  modified  Bessel  functions  and  S is  the  noleculer 
0 \ 
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speed  ratio 

s = 'HP  m 

It  will  be  observed  that  the  drag  force  for  the  non -uniform  case 
differs  from  that  for  the  uniform  eeae  ty  terms  proportional,  respectively , 
to  the  heat  flu*  in  the  X -direction  and  to  the  deviation  of  the  normal 
stresses  from  the  hydrostatic.  The  lift  force,  which  is  zero  for  the  uni- 
form c&uo,  uuu inuio  bci s s proportional  to  the  shear  stress  and  the  heat  flux 
in  the  y direction. 


Eqs, 
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3.1  and  3.2  say  be  rewritten  in  the  form 
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The  indicated  partial  lift  and  drag  coefficients,  , Cpq*  » ate,, 

correspond  to  the  contribution  to  the  total  lift  or  drag  arising  from  the 
non-uniformity  denoted  in  the  subscript.  These  results  are  presented 
graphically  in  Flga.  2 and  3 for  the  case  of  diffuse  reflection.  It  should 
also  be  mentioned  that  the  results  of  Eqs.  3.1  and  3.2  agree  with  those  of 
Ref.  2 for  the  uniform  case. 


(3.4) 

(3.5) 


Inspection  of  Figs.  2 and  3 reveals  that  the  partial  lift  and 
drag  coefficients  corresponding  to  the  noa^uaifora  terms  are  all  small  com- 
pared to  the  coefficient  for  the  uniform  case  C^.  unless  a gas  flow  is 
veiy  non-uniform,  i.e.,  if  the  flow  itself  is  in  the  continuum  or  slip 
flow  region,  -Use  quantities  pxx  , p^y  and  V^y  are  snail  compared  to  p . 

It  follows  that  the  contribution  to  the  total  force  on  the  cylinder  arising 
from  these  viscous  stress  bdisB  wrix  1 be  small  compered  to  the  force  due  to 
the  uniTorm  flow  (of  the  order  of  e few  per  cent  or  Iocs) . Thsss  quantities 
might  thus  be  neglected  in  connection  with  use  of  the  cylinder  aB  a probe. 


The  relative  importance  of  the  forces  arising  from  the  heat  flu* 
terms,  however,  depends  on  the  magnitude  of  the  ratio  of  or  C}y  to  pU  . 
At  high  flow  velocities  these  forces  will  thus  be  negligible  compared  to 
the  force  for  the  unifora  ceae.  At  very  low  velocities  they  will  not.  In 
a stagnation  region  such  forces  constitute  the  entire  force  on  the  cylinder. 
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4.0  CCKCLE5ICNS 

4 = 1 In  the  free  molecule  flow  of  a ncn-unifora  gas  (i.e.,  one 

in  which  heat  flux  or  viecous  "trasses  are  present)  past  a cylinder, 
both  the  lift  and  drag  forces  will  be  affected  by  the  non-unif ornity . 

4.2  Force  contributions  due  to  stress  and  heat  flux  terms  are, 

in  general,  small  compered  to  the  force  for  the  uniform  case,  except 
In  regions  of  low  flow  velocity  where  the  forces  arising  from  tho  heat 
flux  terns  become  1 aperient. 


5.0  mtmciss 

1)  H.  S.  Tsien 

2)  J.  R.  Stalder  and 

V.  J . zurick 

3)  J:  Rr  Stalder, 


G. 

Goodwin  and 

■i 

k*  • 

0.  Ci">>sg6r 

1 * B * 

Bstenfcanij  and 

15. 

D.  Kano 

* W>1 

6)  P.  Epstein 
7}  J.  C.  Maxwell 

8)  H.  Grad 

9)  L.  B.  I<oeb 
10)  M.  Beineman 

U)  H.  Ashley 

12)  C.  S.  Wang  Chang 


■Super*  erodynaaics.  Mechanics  of  Rarefied 
Gases*,  Jour.  Aero.  Sci.,  Vcl.  13,  Dec. 


“Theoretical  Aercx^f..iaic  Characteristics  of 
Bodies  in  a Free  Molecule  Flow  Field®,  SAGA 

TN-2423,  1951 

SA  Cceperiaon  of  Theory  aud  Experiment  for 
High  Spaed  Free  Kolweule  Flow®,  KACA  TN-2244, 
1950 

"A  Torsion  Balance  for  Measuring  Forces  in 
Low  Density  Gab  Flows®,  Jovr.  Applied  Physics, 
Vol.  20,  June  1949,  pp  608-610 

®Zur  Theorie  der  Radiometerkrafte" , 2.  Pbysik, 
Vol.  27,  1924,  p.  1 

■Zur  Theorie  des  Radiometers" , Z.  Phvaik, 

Vol.  54,  1929,  p.  537 

“On  Stresses  in  Rarefied  Gases  Arising  fro* 
Inequalities  of  Temperature" , Trans.  Royal 
Soc  = , Lender.,  Vol.  170,  1879,  p.  231 

"On  the  Einetie  Theory  of  Rarefied  Gases*, 
Coca,  on  Pure  and  Applied  Math.,  Vol,  2, 

1949,  pp  331-407 

■Kinetic  Theory  of  Gases",  McGraw-Hill,  1934 

"Theory  of  Drag  is  Highly  Rarefied  Gases", 
Conns • on  Pure  and  Applied  isath,,  Vol.  1, 

1949,  pp  259-273 

"Applications  of  the  Theory  of  Free  Molseuls 
Flow  to  Aeronautics",  Jour.  Aero.  Sci., 

Vol.  16,  1949,  pp  95-104 

"Transport  Phenooens  in  Very  Dilute  Gases3, 

RCPD  7924-C17H-3-5,  1950 


DRAG  COEFFICIENT 
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